Signatures of adaptive evolution in platyrrhine primate genomes.

The platyrrhine family Cebidae (capuchin and squirrel monkeys) exhibit among the largest primate encephalization quotients. Each cebid lineage is also characterized by notable lineage-specific traits, with capuchins showing striking similarities to Hominidae such as high sensorimotor intelligence with tool use, advanced cognitive abilities, and behavioral flexibility. Here, we take a comparative genomics approach, performing genome-wide tests for positive selection across five cebid branches, to gain insight into major periods of cebid adaptive evolution. We uncover candidate targets of selection across cebid evolutionary history that may underlie the emergence of lineage-specific traits. Our analyses highlight shifting and sustained selective pressures on genes related to brain development, longevity, reproduction, and morphology, including evidence for cumulative and diversifying neurobiological adaptations across cebid evolution. In addition to generating a high-quality reference genome assembly for robust capuchins, our results lend to a better understanding of the adaptive diversification of this distinctive primate clade.

Sensitivity of the polyDetect computational pipeline for phylogenetic analyses.

Alu elements are powerful phylogenetic markers. The combination of a recently-developed computational pipeline, polyDetect, with high copy number Alu insertions has previously been utilized to help resolve the Papio baboon phylogeny with high statistical support. Here, the polyDetect method was applied to the highly contentious Cebidae phylogeny within New World monkeys (NWM). The polyDetect method relies on conserved homology/identity of short read sequence data among the species being compared to accurately map predicted shared Alu insertions to each unique flanking sequence. The results of this comprehensive assessment indicate that there were insufficient sequence homology/identity stretches in non-repeated DNA sequences among the four Cebidae genera analyzed in this study to make this strategy phylogenetically viable. The ~20 million years of evolutionary divergence of the Cebidae genera has resulted in random sequence decay within the short read data, obscuring potentially orthologous elements in the species tested. These analyses suggest that the polyDetect pipeline is best suited to resolving phylogenies of more recently diverged lineages when high-quality assembled genomes are not available for the taxa of interest.

Functional characterization of enhancer evolution in the primate lineage.

BACKGROUND: Enhancers play an important role in morphological evolution and speciation by controlling the spatiotemporal expression of genes. Previous efforts to understand the evolution of enhancers in primates have typically studied many enhancers at low resolution, or single enhancers at high resolution. Although comparative genomic studies reveal large-scale turnover of enhancers, a specific understanding of the molecular steps by which mammalian or primate enhancers evolve remains elusive. RESULTS: We identified candidate hominoid-specific liver enhancers from H3K27ac ChIP-seq data. After locating orthologs in 11 primates spanning around 40 million years, we synthesized all orthologs as well as computational reconstructions of 9 ancestral sequences for 348 active tiles of 233 putative enhancers. We concurrently tested all sequences for regulatory activity with STARR-seq in HepG2 cells. We observe groups of enhancer tiles with coherent trajectories, most of which can be potentially explained by a single gain or loss-of-activity event per tile. We quantify the correlation between the number of mutations along a branch and the magnitude of change in functional activity. Finally, we identify 84 mutations that correlate with functional changes; these are enriched for cytosine deamination events within CpGs. CONCLUSIONS: We characterized the evolutionary-functional trajectories of hundreds of liver enhancers throughout the primate phylogeny. We observe subsets of regulatory sequences that appear to have gained or lost activity. We use these data to quantify the relationship between sequence and functional divergence, and to identify CpG deamination as a potentially important force in driving changes in enhancer activity during primate evolution.

Alu elements and the phylogeny of capuchin (Cebus and Sapajus) monkeys.

Three families of New World monkeys, the Pitheciidae, Atelidae, and Cebidae, are currently recognized. The monophyly of the Cebidae is supported unequivocally by the presence of ten unique Alu elements, which are absent from the other two families. In this paper, the five genomic regions containing these Alu elements were sequenced in specimens representing nine capuchin (Cebus, Sapajus) species in order to identify mutations that may help elucidate the taxonomy and phylogenetic relationships of the cebids. The results confirmed the presence of previously described Alu elements in the capuchins. An Alu insertion present in the Cebidae2 genomic region belonging to the AluSc subfamily was amplified and sequenced only in Sapajus. No amplified or unspecific product was obtained for all other species studied here. An AluSc insertion present in the CeSa1 region was found only in Cebus, Sapajus, and Saimiri. Cebidae4 was characterized by two insertions, an AluSz6 shared by all cebids, and a complete SINE (AluSx3) found only in the capuchins (Cebus and Sapajus). The genomic region Cebidae5 revealed two insertion events, one of the AluSx subfamily, which was shared by all cebids, and another (AluSc8), that was unique to Cebus, offering a straightforward criterion for the differentiation of the two genera, Cebus and Sapajus. The Cebidae6 region showed four distinct insertion events: a 52-bp simple repeat ((TATG) n), two very ancient repeats (MIRc) and a TcMar-Tigger shared by all New World monkeys studied so far, and an Alu insertion of the AluSx subfamily present exclusively in the cebids. The phylogenetic tree confirmed the division of the capuchins into two genera, Cebus and Sapajus, and suggested the southern species Sapajus nigritus robustus and S. cay as the earliest and second earliest offshoots in this genus, respectively. This supports a southern origin for the Sapajus radiation.

From trajectories to averages: an improved description of the heterogeneity of substitution rates along lineages.

The accuracy and precision of species divergence date estimation from molecular data strongly depend on the models describing the variation of substitution rates along a phylogeny. These models generally assume that rates randomly fluctuate along branches from one node to the next. However, for mathematical convenience, the stochasticity of such a process is ignored when translating these rate trajectories into branch lengths. This study addresses this shortcoming. A new approach is described that explicitly considers the average substitution rates along branches as random quantities, resulting in a more realistic description of the variations of evolutionary rates along lineages. The proposed method provides more precise estimates of the rate autocorrelation parameter as well as divergence times. Also, simulation results indicate that ignoring the stochastic variation of rates along edges can lead to significant overestimation of specific node ages. Altogether, the new approach introduced in this study is a step forward to designing biologically relevant models of rate evolution that are well suited to data sets with dense taxon sampling which are likely to present rate autocorrelation. The computer programme PhyTime, part of the PhyML package and implementing the new approach, is available from http://code.google.com/p/phyml (last accessed 1 August 2012).

Chromosome evolution in new world monkeys (Platyrrhini).

During the last decades, New World monkey (NWM, Platyrrhini, Anthropoideae) comparative cytogenetics has shed light on many fundamental aspects of genome organisation and evolution in this fascinating, but also highly endangered group of neotropical primates. In this review, we first provide an overview about the evolutionary origin of the inferred ancestral NWM karyotype of 2n = 54 chromosomes and about the lineage-specific chromosome rearrangements resulting in the highly divergent karyotypes of extant NWM species, ranging from 2n = 16 in a titi monkey to 2n = 62 in a woolly monkey. Next, we discuss the available data on the chromosome phylogeny of NWM in the context of recent molecular phylogenetic analyses. In the last part, we highlight some recent research on the molecular mechanisms responsible for the large-scale evolutionary genomic changes in platyrrhine monkeys.

Retropositional events consolidate the branching order among New World monkey genera.

Due to contradicting relationships obtained from various morphological and genetic studies, phylogenetic relationships among New World monkey genera are highly disputed. In the present study, we analyzed the presence/absence pattern of 128 SINE integrations in all New World monkey genera. Among them, 70 were specific for only a single genus, whereas another 18 were present in all New World monkey genera. The 40 remaining insertions were informative to elucidate phylogenetic relationships among genera. Several of them confirmed the monophyly of the three families Cebidae, Atelidae and Pitheciidae as well as of the subfamily Callithrichinae. Further markers provided evidence for a sister grouping of Cebidae and Atelidae to the exclusion of Pitheciidae as well as for relationships among genera belonging to Callithrichinae and Atelidae. Although a close affiliation of Saimiri, Aotus and Cebus to Callithrichinae was shown, the relationships among the three genera remained unresolved due to three contradicting insertions.

Evolution of the SEC1 gene in New World monkey lineages (Primates, Platyrrhini).

The structure and evolution of the SEC1 gene were examined for the first time in New World primates of the genera Alouatta, Aotus, Ateles, Brachyteles, Callicebus, Callithrix, Cebus, Chiropotes, Lagothrix, Leontopithecus, Pithecia, Saguinus, and Saimiri. This gene has a high CG content (63.8%) and an estimated heterogeneous size ranging from 795 (Callithrix) to 1041 bp (Pithecia), due to numerous indel events. Similar to other fucosyltransferases, three conserved regions are shared by these primates, except for the callitrichines, Aotus and Pithecia, in which indel events resulted in premature stop codons that are related to the production of a supposedly non-functional protein. Phylogenetic analysis of the SEC1 gene, transition/transversion rates, and nucleotide sequence alignment support the hypothesis that primate SEC1 evolved by divergent evolution, and that the lack of activity in some lineages occurred independently at least twice in New World primates, once in the Aotus-Cebus-Callitrichinae group and again in Pithecia. Likelihood-based inference of ancestral states for the activity of SEC1 leads us to suppose that inactivation of SEC1 in the Callitrichinae was a result of a more complex series of events than in Pithecia.

Exploring evolution in Ceboidea (Platyrrhini, primates) by Williams-Beuren probe (HSA 7q11.23) chromosome mapping.

The ancestral platyrrhine karyotype was characterised by a syntenic association of human 5 and a small segment of human 7 orthologues. This large syntenic association has undergone numerous rearrangements in various phylogenetic lines. We used a locus-specific molecular cytogenetic approach to study the chromosomal evolution of the human 7q11.23 orthologous sequences (William-Beuren syndrome, WS) in various Ceboidea (Platyrrhini) species. The fluorescent in situ hybridisation of the WS probe revealed a two-way pattern of chromosomal organisation that suggests various evolutionary scenarios. The first pattern (seen in Callimico and Saimiri) includes a fairly simple disruption of the 7/5 syntenic association by a chromosome fission. The second pattern (seen in Atelinae, Alouattinae and in Callicebus) is characterised by an increasing complexity in the 7/5 association as a consequence of a series of inversions and translocations resulting in different syntenic associations. These data support recent proposals for phylogenomic groupings of New World monkeys. The study also illustrates how single-locus probe hybridisations can reveal intrachromosomal rearrangements.

Variation in maxillary sinus anatomy among platyrrhine monkeys.

Variations in the maxillary sinus anatomy of extant and fossil catarrhine primates have been extensively examined using computed tomography (CT), and have potential utility for phylogenetic analyses. This approach has also been used to demonstrate its anatomical variation in eight of the 16 extant genera of platyrrhines and the absence of the sinus in Saimiri and Cacajao. We used this approach to evaluate the three-dimensional anatomy of the maxillary sinus in all extant platyrrhine genera, and here argue the phylogenic implications of this variation. This study confirms, for the most part, previous CT studies and augments them with the six genera not studied previously: Ateles, Lagothrix, Callithrix, Cebuella, Pithecia and Chiropotes. The entire maxilla is pneumatized by the sinus in the atelines, Cebus, and Callicebus, whereas the sinus pneumatizes only the medial part of the maxilla in the callitrichines and Aotus. Pithecia has a unique conformation in which the maxillary sinus and the expanded inferior meatus pneumatize the posteromedial and anterolateral parts of the entire maxilla, respectively. Chiropotes has no sinus, and the inferior meatus possibly expands into the area between the middle meatus and medial surface of the maxilla to disturb sinus formation, as in the case of its close relative Cacajao. Finally, we argue that the sinus that pneumatizes the entire maxilla is a primitive feature in extant platyrrhines and was probably shared by the last common ancestor of the anthropoids.

The phylogenetic history of New World monkey beta globin reveals a platyrrhine beta to delta gene conversion in the atelid ancestry.

Orthologues of the beta globin gene locus from 10 New World monkey species were sequenced and aligned against available beta and delta globin sequences from rabbit and other primates. Where needed, additional primate sequencing was performed. Phylogenetic analysis identified a beta to delta conversion in the stem of the Anthropoidea, stretching from the 3' part of the proximal promotor to the 5' start of intron 2, consistent with earlier findings. No further conversion appeared to have occurred in the descent of the catarrhines. Within the New World monkey lineage that led to spider monkey and other atelids, another shorter gene conversion was found, spanning adjacent parts of exon 1 and intron 1. The analysis also confirmed that galago beta had replaced galago delta, that an earlier loriform-specific gene conversion extended over intron 2, and that gene conversion throughout the main gene conversion region occurred in the tarsiiform lineage. Platyrrhine phylogenetic relationships were investigated with beta sequences restricted to those that were not involved in gene conversions. This phylogeny generally agreed with results from other nuclear genes. The one exception was that the beta sequences did not place the callitrichine clade within the Cebidae but weakly joined the callitrichine and atelid clades.

Multi-directional chromosome painting maps homologies between species belonging to three genera of New World monkeys and humans.

We mapped chromosomal homologies in two species of Chiropotes (Pitheciini, Saki Monkeys) and one species of Aotus (Aotinae, Owl Monkey) by multi-directional chromosome painting. Human chromosome probes were hybridized to Chiropotes utahicki, C. israelita and Aotus nancymae metaphases. Wooly Monkey chromosome paints were also hybridized to Owl Monkey metaphases. We established Owl Monkey chromosome paint probes by flow sorting and reciprocally hybridized them to human chromosomes. The karyotypes of the Bearded Saki Monkeys studied here are close to the hypothesized ancestral platyrrhine karytoype, while that of the Owl Monkey appears to be highly derived. The A. nancymae karyotype is highly shuffled and only three human syntenic groups were found conserved coexisting with 17 derived human homologous associations. A minimum of 14 fissions and 13 fusions would be required to derive the A. nancymae karyotype from that of the ancestral New World primate karyotype. An inversion between homologs to segments of human 10 and 16 suggests a link between Callicebus and Chiropotes, while the syntenic association of 10/11 found in Aotus and Callicebus suggests a link between these two genera. Future molecular cytogenetic work will be needed to determine whether these rearrangements represent synapomorphic chromosomal traits.

Atelinae phylogenetic relationships: the trichotomy revived?

This research examines phylogenetic relationships between members of the Atelinae subfamily (Alouatta, Ateles, Brachyteles, and Lagothrix), based on analysis of three genetic regions. Two loci, cytochrome c oxidase subunit II (COII) and the hypervariable I portion of the control region, are part of the mitochondrial genome. The other is a single-copy nuclear gene, Aldolase A Intron V. Analysis of these genetic regions provides support for tribe Alouattini containing the Alouatta species, while tribe Atelini contains the other three genera. However, these three genetic regions produce conflicting results for relationships among tribe Atelini members. Previous genetic studies supported grouping Brachyteles with Lagothrix, leaving Ateles in a separate subclade. The present data sets vary based on the genetic region analyzed and method of analysis suggesting all possible cladistic relationships. These results are more consistent with investigations of morphology and behavior among these primates. The primary cause of discrepancy between this study and previous genetic studies is postulated to reside in increased sampling in the present study of genetic variation among members of the Atelinae, specifically Ateles. The present study utilized samples of Ateles from all postulated species for this genetically variable primate, while previous studies used only one or two species of Ateles. This paper demonstrates that shifting relationships are produced when different species of Ateles are used to reconstruct phylogenies. This research concludes that a trichotomy should still be supported between members of tribe Atelini until further analyses, which include additional Atelinae haplotypes are conducted.

Ancestral population sizes and species divergence times in the primate lineage on the basis of intron and BAC end sequences.

The effective sizes of ancestral populations and species divergence times of six primate species (humans, chimpanzees, gorillas, orangutans, and representatives of Old World monkeys and New World monkeys) are estimated by applying the two-species maximum likelihood (ML) method to intron sequences of 20 different loci. Examination of rate heterogeneity of nucleotide substitutions and intragenic recombination identifies five outrageous loci (ODC1, GHR, HBE, INS, and HBG). The estimated ancestral polymorphism ranges from 0.21 to 0.96% at major divergences in primate evolution. One exceptionally low polymorphism occurs when African and Asian apes diverged. However, taking into consideration the possible short generation times in primate ancestors, it is concluded that the ancestral population size in the primate lineage was no smaller than that of extant humans. Furthermore, under the assumption of 6 million years (myr) divergence between humans and chimpanzees, the divergence time of humans from gorillas, orangutans. Old World monkeys, and New World monkeys is estimated as 7.2, 18, 34, and 65 myr ago, respectively, which are generally older than traditional estimates. Beside the intron sequences, three other data sets of orthologous sequences are used between the human and the chimpanzee comparison. The ML application to these data sets including 58,156 random BAC end sequences (BES) shows that the nucleotide substitution rate is as low as 0.6-0.8 x 10(-9) per site per year and the extent of ancestral polymorphism is 0.33-0.51%. With such a low substitution rate and short generation time, the relatively high extent of polymorphism suggests a fairly large effective population size in the ancestral lineage common to humans and chimpanzees.

Species-specific transformation of T cells by HV(MNE).

HV(MNE) is an Epstein-Barr virus (EBV)-like lymphocryptovirus (LCV) originally isolated from a Macaca nemestrina with CD8(+) T cell mycosis fungoides/cutaneous T cell lymphoma (Blood 98 (2001), 2193). HV(MNE) transforms rabbit T cells in vitro and causes T cell lymphoma in New Zealand white rabbits. Here we demonstrate that HV(MNE) also immortalizes T cells from mustached tamarins but not those from owl monkeys, common marmosets, squirrel monkeys, black-capped capuchins, and humans. Cytogenetic and FACS analysis revealed the true origin and T cell lineage of the transformed tamarin T cell lines. Tamarin T cells contained HV(MNE) DNA sequence and displayed a decreased requirement for the IL-2 cytokine for growth. Thus, this EBV-like virus from M. nemestrina differs from the other EBV-like viruses found in nonhuman primates inasmuch as it appears to preferentially transform T cells.

The lowest diploid number (2n = 16) yet found in any primate: Callicebus lugens (Humboldt, 1811).

Morphologic, molecular and karyologic analyses of Callicebus lugens (Humboldt, 1811) of known geographic origin supported the proposition that this is a valid species. Morphologic and morphometric analyses showed evident differences between C. lugens and two other related taxa of the same group (Callicebus purinus and Callicebus torquatus). Cytochrome b DNA analyses (maximum parsimony, neighbour joining and maximum likelihood) were congruent in showing a strong association between C. lugens and Callicebus sp. of the torquatus group in one branch and a sister branch further divided into two clades: one with species of the personatus group and another, with species of the moloch group. Karyotypic analysis showed that C. lugens has the lowest diploid chromosome number of the primate order (2n = 16). Comparisons with other congeneric species clearly supported the proposition that C. lugens is karyotypically similar to others of the torquatus group.

RFLP analysis of mtDNA from six platyrrhine genera: phylogenetic inferences.

This study investigates the phylogenetic relationships of 10 species of platyrrhine primates using RFLP analysis of mtDNA. Three restriction enzymes were used to determine the restriction site haplotypes for a total of 276 individuals. Phylogenetic analysis using maximum parsimony was employed to construct phylogenetic trees. We found close phylogenetic relationships between Alouatta, Lagothrix and Ateles. We also found a close relationship between Cebus and Aotus, with Saimiri clustering with the atelines. Haplotype diversity was found in four of the species studied, in Cebus albifrons, Saimiri sciureus, Lagothrix lagotricha and Ateles fusciceps. These data provide additional information concerning the phylogenetic relationships between these platyrrhine genera and species.

New World monkey phylogeny based on X-linked G6PD DNA sequences.

The Platyrrhini, or New World monkeys, are an infraorder of Primates comprised of 16 genera. Molecular phylogenetic analyses have consistently sorted these genera into three groups: the Pitheciidae (e.g., saki and titi monkeys), Atelidae (e.g., spider and howler monkeys), and Cebidae (e.g., night monkeys, squirrel monkeys, and tamarins). No consensus has emerged on the relationships among the three groups or within the Cebidae. Here, approximately 0.8 kb of newly generated intronic DNA sequence data from the X-linked glucose-6-phosphate dehydrogenase (G6PD) locus have been collected from nine New World monkey taxa to examine these relationships. These data are added to 1.3 kb of previously generated G6PD intronic DNA sequence data [Mol. Phylogenet. Evol. 11 (1999) 459]. Using distance and parsimony-based techniques, G6PD sequences provide support for an initial bifurcation between the Pitheciidae and the remaining platyrrhines, linking Atelidae and Cebidae as sister taxa. Bayesian methods provided a conflicting phylogeny with Atelidae as outgroup. Within the Cebidae, a sister relation between Aotus and the Cebus/Saimiri clade is favored by parsimony analysis, but not by other analyses. Potential reasons for the difficulty in resolving family level New World monkey phylogenetics are discussed.

Molecular cladistic markers in New World monkey phylogeny (Platyrrhini, Primates).

Transpositions of primate-specific Alu elements were applied as molecular cladistic markers in a phylogenetic analysis of South American primates. Seventy-four human and platyrrhine loci containing intronic Alu elements were PCR screened in various New World monkeys and the human outgroup to detect the presence of orthologous retrotransposons informative of New World monkey phylogeny. Six loci revealed size polymorphism in the amplification pattern, indicating a shared derived character state due to the presence of orthologous Alu elements confirmed by subsequent sequencing. Three markers corroborate (1) New World monkey monophyly and one marker supports each of the following callitrichine relationships: (2) Callithrix and Cebuella are more closely related to each other than to any other callitrichine, (3) the callitrichines form a monophyletic clade including Callimico, and (4) the next living relatives to the callitrichines are Cebus, Saimiri, and Aotus.

A new subclass of the luteinizing hormone/chorionic gonadotropin receptor lacking exon 10 messenger RNA in the New World monkey (Platyrrhini) lineage.

The luteinizing hormone receptor (LHR) plays an essential role as a mediator of LH and CG action during embryonic sexual differentiation and in gametogenesis. In a hypogonadal male patient, we recently demonstrated that a genomic deletion of exon 10, located in the hinge region of the extracellular domain, results in discrimination of LH and hCG action. In the common marmoset (Calltithrix jacchus), exon 10 of the LHR is naturally missing at the mRNA level. In order to investigate whether this is an isolated species-specific phenomenon, we performed a phylogenetic screening, searching for the presence of LHR exon 10 mRNA in a number of primate species representative for the major lineages of primate evolution. The expressed LHR region encompassing exon 10 was amplified from testicular tissue by RT-PCR, cloned, and sequenced. In addition, we performed Southern blot analysis of the LHR of selected New World and Old World primates. The results revealed that exon 10 mRNA is lacking in the complete New World monkey (Platyrrhini) lineage but is present in both more primitive and more advanced primates. However, exon 10 seems to be present at the genomic level, arguing for a splicing failure possibly due to a genomic mutation or the lack of appropriate splicing factors. Considering that, in the human, LH is far less active than hCG on the LHR lacking exon 10, we addressed the question whether the existence of such a receptor has any consequences on the dual hormone LH/CG system present in Platyrrhini. Using primers specific for the known marmoset CG beta cDNA, we amplified the CG beta subunit cDNA from male common marmoset pituitaries by RT-PCR, while LH beta could not be amplified, suggesting a possible physiological role of pituitary CG in this species. In conclusion, we demonstrated for the first time that the LH mRNA without exon10 is the natural wild-type LHR in the Platyrrhini lineage. We propose that this LHR represents a new subclass of receptors that should be named LHR type II. In addition, the high expression of CG beta in the marmoset pituitary suggests a physiological role of CG in the reproductive function of these primates beyond pregnancy.